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We study various properties of an ultracold two-dimensional (2D) Bose gas that are beyond a 
mean-field description. We first derive the effective interaction for such a system as realized in 
current experiments, which requires the use of an energy dependent T- matrix. Using this result, 
we then solve the mean-field equation of state of the modified Popov theory, and compare it with 
the usual Hartree-Fock theory. We show that even though the former theory does not suffer from 
infrared divergences in both the normal and superfluid phases, there is an unphysical density discon- 
tinuity close to the Berezinskii-Kosterlitz-Thouless transition. We then improve upon the mean-field 
description by using a renormalization group approach and show how the density discontinuity is 
resolved. The flow equations in two dimensions, in particular, of the symmetry-broken phase, al- 
ready contain some unique features pertinent to the 2D XY model, even though vortices have not 
been included explicitly. We also compute various many-body correlators, and show that correlation 
effects beyond the Hartree-Fock theory are important already in the normal phase as criticality is 
approached. We finally extend our results to the inhomogeneous case of a trapped Bose gas using the 
local-density approximation and show that close to criticality, the renormalization group approach 
is required for the accurate determination of the density profile. 



INTRODUCTION 



Low-dimensional systems play a unique role in the 
study of many-body effects. For example, the en- 
hanced importance of thermal fluctuations prevents a 
two-dimensional (2D) system with a continuous symme- 
try to undergo spontaneous symmetry breaking at any 
nonzero temperature, thus preventing the presence of 
true long-range order. This property is elucidated in the 
Mermin-Wagner-Hohenberg theorem [1, 2]. Nonetheless, 
the system still exhibits interesting properties. In par- 
ticular, the 2D XY model undergoes a special type of 
phase transition into a state which is characterized by 
only algebraic long-range order instead. The underly- 
ing mechanism that drives the phase transition, known 
as the Berezinskii-Kosterlitz-Thouless (BKT) transition 
[3, 4], is the unbinding of vortex-antivortex pairs. Due 
to its topological nature, such a phase transition is dif- 
ficult to incorporate into the standard Ginzburg-Landau 
theory with a local order parameter. For the ultracold 
2D Bose gas, the absence of Bose-Einstein condensation 
(BEC) requires the concept of a quasi-condensate to un- 
derstand the existence of algebraic long-range order and 
superfluidity [5-9]. Experiments in the field of ultracold 
atomic gases have recently reached this interesting 2D 
regime to allow for the direct observation of this phe- 
nomenon in a highly controllable environment [10-14]. 
The observation of dislocations in the interference pat- 
tern of two condensates [10, 12, 15, 16] and the studies 
of the coherence properties [11, 17], for example, have 
all agreed with the universal predictions of BKT theory. 
More experiments and numerical simulations have come 
to address also various nonuniversal properties specific to 
the atomic gases under consideration [13, 18-21]. 

To describe the low-dimensional Bose gas at nonzero 
temperatures, the usual approach of the Bogoliubov the- 



ory is plagued with infrared divergences. However, these 
divergences can be shown to occur due to a spurious con- 
tribution from the condensate phase fluctuations in the 
Bogoliubov approach and by removing these contribu- 
tions we can arrive at a modified Popov theory, which 
is valid for any dimension and at all temperatures [8]. 
In this paper, we first study this modified Popov theory 
for the ultracold 2D Bose gas and show that it contains 
a density discontinuity above the BKT transition, close 
to the point where the quasicondensate density becomes 
nonzero. To improve upon the mean-field description, we 
next employ a renormalization group (RG) approach to 
take into account the quantum and thermal fluctuations 
more accurately, in particular in the normal phase, when 
the modified Popov theory reduces to Hartree-Fock the- 
ory. The RG theory developed for the ultracold three- 
dimensional (3D) Bose gas was shown to be quantita- 
tively successful in addressing effects beyond mean-field 
[22]. Here, we derive the analogous 2D RG flow equa- 
tions and find the surprising result that they show fea- 
tures which are very different from the 3D RG theory. We 
first show how various characteristics of a quasiconden- 
sate are manifested in this framework. We then interpret 
these unique features as precursors of the BKT physics, 
even though we have not taken topological defects into 
account explicitly. This is because the RG theory is de- 
rived from the full atomic quantum field ?/;(x, r) and not 
from its phase alone [23]. With the RG approach, the 
density correction to the Hartree-Fock theory indeed re- 
solves the unphysical density discontinuity in the mean- 
field description. Furthermore, it agrees with the equa- 
tion of state of the modified Popov theory already above 
the critical temperature, which shows that the latter has 
correctly included correlation effects beyond the Hartree- 
Fock description even in the normal phase. 

The paper is organized as follows. In Sec. II, we 
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discuss the exact form of the T-matrix for ultracold 2D 
gases as reahzed under current experimental conditions, 
since it determines the effective interaction of the 2D 
Bose gas. Next, in Sec. Ill, we present the mean- field 
results from the modified Popov theory, and point out 
that the Hartree-Fock theory becomes unstable already 
above the BKT transition, even though a solution to 
the Hartree-Fock equation of state exists at all tempera- 
tures. The discontinuity in the density close to the BKT 
transition that results from the modified Popov theory 
is deemed unsatisfactory, and in Sec. IV, we present the 
required RG theory to improve upon this result. We 
discuss various interesting features of the flow equations 
in two dimensions and compute various nonuniversal 
quantities of interest within this approach. While the 
RG approach resolves the density discontinuity of the 
mean-field equation of state, it remains incapable of 
capturing all the critical properties known from the 
BKT physics. We also compute various many-body 
correlators, and show how correlation effects beyond a 
mean-field picture show up in the normal phase close to 
criticality. In Sec. V, we extend our results to the case 
of a trapped Bose gas and compare them with the ideal 
gas. We end with some concluding remarks in Sec. VI. 



(2 7r)-^/^( 2^n! l)-^^^Hn{z/V2l) exp(-zV4/^), where / = 
yjhjmujz is the harmonic length in the axial direction 
and Hn{z) is the Hermite polynomial. 

The Lippmann-Schwinger equation for the two-body 
T-matrix is given by 



f(£;+) = y + y 



(2) 



where = E -\- iO and the effective interaction in the 
2D regime is given by the matrix element with respect to 
the axial ground state, i.e.. 



Too{E) = (k, n = 0|f (^)|k^ n = 0). 



(3) 



By inserting the completeness relation 
J]]^ J (ik|k, n)(k, n|/(27r)^ = 1, the operator equa- 
tion in Eq. (2) can be written as 



T{E) 



-J OO 

7-E 



n=0 



{2TTy E-Eo{k, n)' 



(4) 



where T(E) is the matrix element with respect to 3D 
plane waves, which is related to the desired quantity 
Too{E) by means of Too (^) = |^o(0)pT(£;). Moreover, 
Eo(k,n) are the eigenvalues of the free Hamiltonian 



II. EFFECTIVE INTERACTION IN THE 2D 
REGIME 

Even though we shall be interested in ultracold 2D 
Bose gases, in experiments with atomic alkali- metal gases 
the realization of such a system is achieved by restrict- 
ing the motion of a trapped 3D gas onto a plane. For 
the quantum degenerate gas to be in the 2D regime, one 
needs to ensure that the motion along the tightly confin- 
ing axial direction is frozen out. This condition is met 
if /c^T, /i <^ hujz, where ksT is the thermal energy, 
is the chemical potential, and ojz is the axial trapping 
frequency. 

In the ultracold limit of the gases under considera- 
tion, the effective interaction is in the first instance deter- 
mined by the three-dimensional two-body T-matrix. We 
therefore begin with considering the full two-body Hamil- 
tonian Ho -\- V in the center-of-mass coordinate frame, 
where the free Hamiltonian ^o is given by 
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m is the atomic mass, and V is the interaction poten- 
tial modeled by a short-ranged delta function of strength 
Vq. The center-of-mass coordinate frame is spanned 
by the r-plane, which is taken to be homogenous, and 
the z-axis in the tightly-confining direction. We de- 
note the eigenstates of by |k, n), which are given by 
the product of 2D plane waves (r|k) = e^^'^ and one- 
dimensional oscillator functions {z\n) = ^n(^)- The lat- 
ter functions are defined in the standard form, 0n(^) = 



^o(k,n) = 2ek 



(5) 



with ek = /i^k^/2m. 

As it stands, Eq. (4) suffers from an ultraviolet di- 
vergence due to the delta function interaction potential, 
which neglects any momentum dependence at high mo- 
mentum. To cure this divergence, we observe that while 
the harmonic oscillator functions with odd n vanish at 
the origin, the functions with even n have the asymp- 
totic behavior 
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for large quantum numbers n, where T{n) denotes the 
Gamma function. It then follows that, for large n, the 
second term on the right-hand side of Eq. (4) takes the 
form 
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(27r)3 E-2ek-2efc, 



(7) 



with kz = v2n//. This form diverges in the ultraviolet 
in exactly the same manner as the two-body T-matrix in 
the homogenous three-dimensional space, i.e., as 
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with q = (k, kz). On physical grounds, this behavior is 
expected since the short-distance behavior of the gas is 
not altered by the harmonic confinement and a standard 
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renormalization procedure can therefore be implemented. 
The latter is achieved by eliminating the bare coupling 
strength Vq in favor of the relevant physical parameter 
T'^^(E :^ 0) = Anh'^a/m for the low-energy scattering, 
where a is the 3D scattering length. After the elimina- 
tion, a well-defined expression for the T-matrix, 
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that no longer contains any divergence is obtained. It is 
shown in the Appendix that the effective interaction can 
be further worked out to yield 



tiD 1.0 

0.6 



0.2 



3 4 
E/2hO), 



FIG. 1: The exact T-matrix element |Too(^)| as a function of 
energy. 
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where Fq c± —0.3508, and for convenience energies are 
measured with respect to the zero-point energy 
of the axial ground state. This result can be compared 
with the previous work of Petrov et al. [24], where a 
different approach has been adopted. Thus, we see that 
the interaction takes a quasi-2D form, even though the 
system is kinematically two dimensional. The T-matrix 
Too{E) interpolates between the form of the interaction 
for the 3D and 2D cases, where the former is a constant 
and the latter has a logarithmic energy dependence. For 
the trapped atomic gases of interest, the condition l/a ^ 
1 is usually satisfied. Hence the effective interaction is 
often well approximated by Toq{E) o:^ 2^/27rafi^ /ml = 
g. Nonetheless, the exact form of the T-matrix with an 
energy dependence will turn out to be important in the 
present work. In Fig. 1, we show |Too(-E')| as a function of 
energy. Besides that the overall scale is set by ^, there are 
zeros at energies 2nhwz corresponding to the harmonic 
oscillator states with even n, due to the divergence of the 
logarithmic function. Note again that we have subtracted 
the zero-point energy. 



III. MEAN-FIELD THEORY 

In this section, we recapitulate the main result ob- 
tained from the modified Popov theory for the ultracold 
2D Bose gas [8] . The equation of state in this mean- field 
theory is given by 



n = no + 



j27V(a^k) + l]-2 



1 , noToo(-2/i) 



2ek + 2/i 



and II = (2n — no) Too (—2/i), where n is the total density, 
no is the (quasi) condensate density, n' is the density fluc- 
tuations around no, hjj\^ = [e^ + 2noToo(— 2/i)ek]^/^ is the 
Bogoliubov quasiparticle dispersion, N{x) = l/(e^^ — 1) 
is the Bose-Einstein distribution function, and (3 = 
l/ksT is the inverse thermal energy. Notice that these 
equations do not suffer from infrared and ultraviolet di- 
vergences, and thus, they are valid for any dimension and 
at all temperatures. This comes about because the con- 
densate phase fluctuations have been treated exactly to 
arrive at this result. A condensate with true long-range 
order is absent in two dimensions at any nonzero temper- 
ature. In fact, as shown in Ref. [8], no should be iden- 
tified with the quasicondensate density at any nonzero 
temperature. Due to its mean-field nature, however, this 
equation of state is incapable to capture the BKT tran- 
sition. Indeed, although the criterion for the BKT tran- 
sition is n^A^j^ = 4, where Ath = \/2i^W' jmkBT is the 
thermal de Broglie wavelength of the atoms in the gas, 
according to Eq. (12), a nontrivial solution exists even if 
the superfluid density n^ obeys n^A^^^ < 4. To circumvent 
this shortcoming, it was shown that by complementing 
the above equation of state with a renormalization group 
analysis on the Sine- Gordon model that takes into ac- 
count vortices explicitly, the canonical criterion for the 
BKT transition is met when noA^^^ = 6.65 [8]. Thus, 
above this critical temperature, the fugacity of the vor- 
tices renormalizes at long wavelengths to a nonzero value, 
which leads to the destruction of superfluidity but not 
immediately of the quasi-condensate. 

The normal state without a quasicondensate is instead 
described by the Hartree-Fock equation of state given by 



where the Hartree-Fock self-energy ftS satisfies 
ftS = 2nToo(-/iS). 



(13) 
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= no-|-n' 
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It is also noted that when the T-matrix is energy indepen- 
dent, i.e., Too(^) — g-> the normal equation of state can 
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FIG. 2: The mean- field equation of state from the modified 
Popov theory. The dashed curve is the Hartree-Fock theory 
that is extended also into the unstable regime until the Monte 
Carlo critical condition ric A^h = 7.81 is met (dotted line). We 
have taken mg/h^ = 0.15, which roughly corresponds to the 
experiment of Ref. [10]. The canonical BKT criterion derived 
from the modified Popov theory is indicated with an asterisk. 



actually be satisfied at all temperatures, since Eq. (13) 
can then be integrated to give 



■ln(l 



-(3{2ng-^i) 



(15) 



For a fixed total density n, Eq. (15) has always a solution 
for the chemical potential //, at all temperatures. Thus 
to ensure the stability of the normal phase, we have to 
examine with the same chemical potential ji whether the 
solution of the modified Popov theory exists. 

To illustrate this discussion, we numerically solve both 
equations of state, as shown in Fig. 2. We observe that 
there is an unphysical discontinuity in the density curve 
due to a mismatch in the total density between the two 
equations of state close to the BKT transition. This is an 
artifact of the mean-field theory which will be resolved 
in the next section. 

To compare, we include in Fig. 2 the result from 
the usual Hartree-Fock mean-field approach. In this ap- 
proach, the Hartree-Fock theory is employed up to a 
chemical potential that corresponds to the critical den- 
sity Uc obtained from the high precision Monte Carlo 
simulations [25]. Numerical simulations yielded a crit- 
ical density for the onset of the BKT transition of 
ncA^j^ ln(380 fi^/m^) = 7.81, and a critical chemical 
potential iic/kBT {mg /7ih^)\n{l?>.2h^ /mg) = 0.22, 
for mg/h? c± 0.1541. Both conditions are shown by the 
dotted lines in Fig. 2. It is noted that the critical chem- 
ical potential obtained within this mean-field approach, 
fJ^c^/kBT c± 0.38, differs from the Monte-Carlo result, 
thus indicating an inconsistency in the Hartree-Fock ap- 
proach. This can be understood within the modified 
Popov theory as an instability of the system toward a 
more correlated phase. On the other hand, despite the 
density discontinuity, the modified Popov theory yields 
a critical density n^^ :^ 7.94 and a critical chemical 




FIG. 3: The RG flow diagram of Vb and /x in the symmet- 
ric phase. The arrows indicate the flow direction towards the 
long- wavelength regime. The dotted lines are ill-defined tra- 
jectories as explained in the text. 



potential fif^/ksT 0.23, which are in excellent agree- 
ment with the Monte Carlo results. This shows that the 
main problem of the modified Popov theory lies in the 
inaccurate treatment of correlations in the region where 
it reduces to Hartree-Fock theory. 



IV. RENORMALIZATION GROUP THEORY 

To improve upon the modified Popov theory, we em- 
ploy here an RG approach [26]. In this approach, we 
systematically integrate out the high-momentum shell 
Ae~* < A: < A, and absorb its contribution into the pa- 
rameters of the theory, which then become dependent on 
the flow parameter t. Here, A is the ultraviolet cutoff 
of the theory that is specified below. As discussed in 
Ref. [22], due to the ultracold limit of the Bose gas un- 
der consideration, the parameters which are important 
in determining the various properties of interest are the 
chemical potential /i = /i(t), and the two-body interac- 
tion strength Vb = Vo{t). 



A. The Flow Equations 

We shall present here the flow equations for these pa- 
rameters in two dimensions, thus extending the results 
obtained in Refs. [22, 27]. In the symmetric phase, for 
/i < 0, the flow equations are given by 
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dt 



2/i 



A2 



A2 
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1 + 27V(eA - /i) 
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(16) 
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FIG. 4: The RG flow diagram in the symmetry-broken 
phase. There is a limiting curve which ends at the attract or 
(0, — eA/2) where all trajectories approaches asymptotically. 
The inset shows the approach to the attractor. 



while in the symmetry-broken phase, for /i > 0, they are 



dt 



A2 



2el + Gfiel + fi^ 



[2N{hujK) + 1] 



A2 



(eA - m)^ 



2h?u 



[2N{fiWK) + 1] 



+ 



(2eA + ^if 



l3N{nj^K)[N{tiL0K) + l] 



(17) 



Here, the chemical potential has been trivially rescaled 
by /i ^ ^e^^ hence the 2/i term in the flow equations. 
The inverse temperature has also been trivially rescaled 
as /5 ^ but the interaction strength does not 

scale trivially in two dimensions. To examine the crit- 
ical properties, however, a different trivial scaling needs 
to be used. This comes about because close to the critical 
regime, where the correlation length and correlation time 
diverge, the time-derivative term in the quantum action 
can be neglected with respect to the kinetic term. This 
is equivalent to taking the large-t or high-temperature 
limit of the flow equations by setting N{x) 1/ [3x. As 
a result, while the trivial scaling of the chemical potential 
remains the same, the coupling strength acquires a trivial 
scaling with exponent 2, i.e., Vb ^ Vbe^*. As expected, 
the trivial scalings in the critical regime then agree with 
the trivial scalings of the classical 2D XY model. 

To solve the flow equations, the correct boundary con- 
ditions have to be provided. The initial value for the 
chemical potential = 0) is nothing but the bare chem- 
ical potential in the original theory. Its initial value can 
either be positive or negative, corresponding to an ini- 
tial phase which is symmetry broken or unbroken, re- 
spectively. Furthermore, the flow equations allow the 
chemical potential to change sign, since the two sets of 
equations are smoothly connected at /i = 0. As we shall 
see, it is in fact a general feature of the solutions to the 
2D flow equations in the symmetry-broken phase that 



FIG. 5: The intercept with the Vb axis as a function of the 
chemical potential fi. 



the chemical potential always flows to a negative value 
for large t, which is not the case in three dimensions. 

To determine the boundary condition for the interac- 
tion strength, we recognize that in vacuum, i.e., by set- 
ting A^(eA — /i) = 0, the flow equation for the interaction 
strength 



dt 



A2 



1/2 



2^2(eA-/i) 



(18) 



is nothing but the differential form of the Lippmann- 
Schwinger equation for the T-matrix at energy 2//. Since 
the T-matrix solution to the Lippmann-Schwinger equa- 
tion, as obtained in Sec. II, entails the summation of all 
ladder diagrams for the scattering process in vacuum, we 
have to ensure that the flow equation in Eq. (18) repro- 
duces the correct long- wavelength result for large t. In 
particular, the initial value V^it = Q) is chosen such that 
for large t the correct form of the T-matrix is recovered 
in the vacuum, i.e., V{){t oo) = Too (2/i). This can be 
satisfied by the following initial condition: 



Fo(t = 0) = 



1 



m l/a+ (1/V^) In {mu^/2hK^) 



for f?lS? jm ^ 2/i. It is important to note that for 
A ^ Ath, this procedure allows us to eliminate the 
ultraviolet cutoff dependence of the theory since the 
interaction strength already attained the two-body 
T-matrix value before entering the thermal regime [22]. 
With these boundary conditions and the ultraviolet 
cutoff chosen to be A ~ O (100 Ath), we now numerically 
integrate the flow equations to obtain the specific 
solution {Vb(t), /i(t)} for different chemical potentials. 



B. Analysis of the Flow Equations 

We first study the two sets of flow equations in Eq. (16) 
and Eq. (17) separately, and allow for the chemical po- 
tential in both equations to take positive and negative 
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FIG. 6: The RG flow diagram of the two sets of flow equations 
in the respective regions. 



values. In the symmetric phase, as seen in Fig. 3, the tra- 
jectories seem to be characteristics of the order- disorder 
phase transition of the classical 3D XY model. The fixed 
point of the flow equations can be found by considering 
their large-t limit 



~dt 
dt 



2Vo 



= 0, 



TT 



ex- II 



27r (eA-/i)2 



(20) 



which yields (Vo*,M*) = {2{)Tiel/A^K^kBT,2eK/l). How- 
ever, it is important to note that for chemical poten- 
tials above a critical value, the resulting trajectories are 
actually ill-deflned because eventually grows to a 
point where the Bose-Einstein distribution diverges, i.e. 
7V(eA — /i) = oo. Thus, there is not really an order- 
disorder phase transition with increasing chemical po- 
tential. More generally, the flrst quadrant with positive 
chemical potential should be regarded as the unphysical 
region of the symmetric phase because the true minimum 
of the action is shifted away from the origin. Similarly, 
in three dimensions the unstable flxed point lying in this 
region does not give the most accurate critical properties 
of the system [22]. 

In the symmetry-broken phase, the flow equations 
present a few rather surprising features, as shown in Fig. 
4. First, all trajectories always flow into the fourth quad- 
rant with a negative chemical potential. This feature can 
be interpreted as the prima facie property of a quasicon- 
densate in two dimensions, where starting with an initial 
"condensate" (/i > 0) at the shortest distance, as the high 
momentum shells are being integrated out, the symme- 
try always gets restored (// < 0) in the infrared regime. 
In other words, the high- momentum fluctuations at long 
wavelength always destroy the coherence present at short 
distances. 

Second, there exists an attr actor in the fourth quad- 
rant towards which all the trajectories flow, without the 
need to flne tune the initial condition. The attractor can 
be found from solving the large-t limit of the flow equa- 
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FIG. 7: The RG equation of state, with mglf? = 0.15. The 
dashed lines are the critical conditions for the BKT transition 
from the Monte Carlo (MC) results. The dotted line is the 
modified Popov theory. 



tions. 
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(^o*,M*) 



[eA(eA 
f 2eAfi ■ 



[eA(eA + 2/i)]2 



= 0, 



0, (21) 



(0, — eA/2). Due to a non- 
analytic behavior of the differential equations near the 
attractor, an expansion of Eq. (21) around the attractor 
(Vo*+^^o,/i* + ^/i) gives 



dSjj. 

~dr 



ddV^ 
dt 



IbA^ksT SVo 



dA^keTeA SVp 
SVo, 



A^ksT 

IGtt Sfi Sttca 



(22) 



which thus does not permit a linearization. The approach 
to the attractor can nevertheless be found by substituting 
the ansatz SVq (x 1/t^. In this manner we flnd that the 
direction of approach is given by 



t 



327r 1 + t 



SVoit), 



(23) 



and S/LL oc 1/t for large t. This is shown in the inset of 
Fig. 4. 

Third, there exists a limiting curve that all trajectories 
approach asymptotically. The intercept of the trajecto- 
ries with the Vq axis for increasing chemical potential is 
shown in Fig. 5. The interesting feature here is that for 
sufficiently large chemical potential the intercept essen- 
tially remains constant. We will come back to this point 
shortly. 

Again, since the flow equations in the symmetry- 
broken phase are valid for /i > 0, the fourth quadrant 
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FIG. 8: Superfluid fraction (solid curve) as a function of tem- 
perature for a density of 1.5 x lO^cm"^, with mg/h^ = 0.15. 
The dotted hne extrapolates the solid curve from the RG cal- 
culation to the lower temperature regime. The critical BKT 
condition is shown by the dashed line. 



should be considered as the unphysical region. How- 
ever, the existence of an attr actor and a limiting curve in 
the solution to the flow equations already suggests fea- 
tures which are unique to the ultracold 2D Bose gas, even 
though we did not take vortices into account explicitly. 
They are in fact related to the fixed line known from the 
BKT theory. 

Finally, let us consider the system of coupled flow equa- 
tions, and use the appropriate ones in the different quad- 
rants. The result is shown in Fig. 6. The main point 
to note here is that a critical chemical potential /ic can 
be identified for a fixed temperature, beyond which the 
resulting flow for large t is insensitive to its initial value. 
This is so because for fi > fic^ the trajectory changes the 
sign of at essentially the same value of Vb, and thus 
the continuing flow into the fourth quadrant governed by 
the symmetric phase flow equations has essentially the 
same initial value. In other words, the long- wavelength 
action takes the same form for all fi > fic As a result, we 
identify this /ic with the critical condition for the BKT 
transition in this approach. It is interesting to note that 
the critical chemical potential obtained in this manner 
agrees within 10% with the Monte Carlo result. 



C. RG Equation of State and Correlation Effects 

Having discussed the properties of the flow equations, 
we next determine various nonuniversal quantities of in- 
terest. With the RG approach, the total density n and 
the superfluid density Us can be computed by integrating 
these quantities along the trajectory, which are expressed 
by the following differential equations 



dn 
~dt 
dn 
~dt 



A! 
A! 

27r 



N{eA - for /i < 0, 



e for /i > 0, 



(24) 



dus 
~dt 



dn 
~dt 



A: 

27r 



eKPN{x) [N{x) + 1] e 



-2t 



(25) 



where x = hjO\ or {e\ — fi) for /i > or /i < 0, respec- 
tively. The initial conditions are n{t = 0) = ^^(t = 0) = 
0. In Fig. 7, the density curve n(t ^ oc) as a function of 
the chemical potential is shown and compared with the 
modified Popov theory. For small chemical potentials, 
the small deviation from the mean-field theory is consis- 
tent with the expectation that at low densities the RG 
correction is unimportant. As the chemical potential ap- 
proaches the BKT critical point, however, the deviation 
from the Hartree-Fock theory becomes substantial. In 
fact, the RG density curve connects smoothly with the 
density curve obtained from the modified Popov theory. 
Thus, the RG approach resolves the artificial disconti- 
nuity observed in the mean- field theory. Furthermore, it 
shows that the equation of state for the quasi-condensate 
provides already a good description above the critical 
temperature. 

Next, by integrating Eq. (25), the superfluid density 
is obtained, which is shown in Fig. 8. We see that the 
anticipated discontinuous jump in the superfluid density 
is absent at ^s^th ~ which shows that not all univer- 
sal features of the BKT transition are incorporated yet. 
We believe that the effects on the superfluid density as- 
sociated with the proliferating vortices can be taken into 
account explicitly by performing an additional renormal- 
ization group analysis on the Sine-Gordon model, as done 
in Ref. [8]. The initial condition for the dielectric con- 
stant K{0) = jSh^ns/m should then be identified with 
the superfluid density obtained here. 

Finally, we compute various many-body correlators 
where enhanced correlation effects are expected to show 
up as criticality is approach. From the modified Popov 
theory, the renormalized density- density correlator is 
given by [8] 



/\21 



(26) 



The reduction in the three-body recombination is given 

by 



-N 



L{T) 



Too (-2m) 



Tm{-2hT,) 



Kf{T) 



(27) 



where is the recombination rate constant in the nor- 
mal phase, and the renormalized three-body correlator is 
given by 



1 



9ngn + 18no(n )^ + Q{n'f] . (28) 



In Fig. 9 and Fig. 10, we see that as the system in the 
normal phase approaches criticality, the mean-field de- 
scription of Hartree-Fock theory for the many-body cor- 
relators breaks down completely. While the latter gives a 
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FIG. 9: Renormalized density-density correlator as a function 
of chemical potential, with mg/h^ = 0.15. The dotted line 
is the prediction from the Hartree-Fock theory, with (n^) = 
2(n)^. The dashed line is the BKT transition from the Monte 
Carlo result. 
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FIG. 10: Reduction of the three-body recombination rate as 
a function of the density at temperature T = 84.4 nK. The 
dashed line is the BKT transition from the Monte Carlo re- 
sult. 



constant value of 1 for the renormalized density-density 
correlator, the enhanced correlation due to the presence 
of a quasicondensate reduces this value to 0.68 at criti- 
cality. Furthermore, the three-body recombination rate 
reduces by a factor of 2.8, compared to the recombination 
rate constant in the normal phase. Thus, both these can 
serve as an observable for beyond Hartree-Fock effects. 



V. TRAPPED BOSE GASES 

In this section, we extend our results to the inhomo- 
geneous case of trapped Bose gases. For the case of a 
trapped ideal Bose gas, the phenomenon of Bose-Einstein 
condensation (BEC) occurs. Within the local-density ap- 
proximation (LDA), BEC takes place when the phase- 
space density in the center of the trap diverges. The 
temperature at which this phenomenon occurs, the BEC 
temperature Tbec^ can be related to the total particle 
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FIG. 11: Density profiles for a quantum degenerate cloud 
of ^^Rb atoms for various theories for N = 7410 at T = 
Tbec- We take a = 5.2nm, a radial trapping frequency of 
uj — y/uJxOJy — 27rA/9.4 x 125 Hz ^ 27r x 34.3 Hz and an axial 
trapping frequency a;^ = 27r x 4000 Hz. With this particle 
number, the phase-space density in the center of the trap 
diverges for the ideal Bose gas. 

number N in the trap by [28] 

where O is the geometric mean of the radial trapping 
frequencies. For the case of a trapped interacting Bose 
gas in two dimensions, the system does not undergo a 
BEC, but a BKT transition if the trapping frequency O 
is sufficiently low [29, 30]. For current experiments of in- 
terest, the harmonic length ^Jh|m^ associated with the 
radial trapping frequencies easily exceeds the de Broglie 
wavelength Ath- Thus, the use of LDA by incorporating 
the effect of radial trapping through the introduction of 
a local chemical potential /i(r) = /i — mcj^r^/2 is readily 
justified. We show in Fig. 11 the density proffie of an 
ideal Bose gas at the critical condition and compare it 
with the case of an interacting Bose gas within the LDA. 
Here, we see the drastic effect of interactions in the 2D 
Bose gas. For these conditions, the RG density profile 
for the interacting gas is only slightly different from the 
Hartree-Fock theory, which is stable in this case. We in- 
clude the density profile of an ideal classical Boltzmann 
gas for comparison. 

As the phase-space density is increased, the Hartree- 
Fock theory is no longer applicable due to the instability 
towards a more correlated state with a quasicondensate, 
as discussed previously. We show the RG density proffies 
in this range of phase-space density in Fig. 12. The two 
density profiles describe two different regimes of the gas, 
namely, one with the phase-space density in the center 
of the trap slightly below the critical phase-space density 
ncA^j^ for the BKT transition, and the other well above 
it. We note that in the latter case, the gas consists of a 
superfiuid core up to a critical radius, and is surrounded 
by an outer normal shell. One quantity of interest is the 
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FIG. 12: Density profiles from RG at T = 0.83 Tbec and T = 
O.GITbec, for N = 10800 and TV = 20000, respectively. In 
this range of phase-space density, the Hartree-Fock theory is 
unstable towards the equation of state of the modified Popov 
theory. The dash-dotted line is the critical condition for the 
BKT transition. 



BKT transition temperature in the trap, relative to the 
ideal gas BEG temperature. By integrating the density 
profile at criticality, we obtain Tbkt — 0.81 Tbec for 
mg/fi^ = 0.15. 

There has been recent interesting experimental and 
numerical work [13, 18-21] being carried out to address 
the nonuniversal quantities of the quasi-2D trapped Bose 
gases. A direct comparison with our theory is rendered 
difficult, because in these works, the thermal excitations 
in the tightly confining direction are non- negligible. 
Throughout our paper, we consider, instead, the strictly 
2D regime, i.e., /i, ksT <C h(jJz- 



We then presented the mean-field results of the modified 
Popov theory. Even though the theory can describe 
both the normal and superfluid states, the density turns 
out to be discontinuous close to the BKT transition, 
which is an artifact of the theory. We improved upon 
the mean-field description by a RG approach. The 
flow equations exhibit interesting features, which re- 
semble many of the unique properties of the 2D XY 
model, even though the effects of vortices have not 
been included explicitly. With the RG approach, the 
density correction to the normal equation of state indeed 
connects smoothly with the quasicondensate equation 
of state in the superfluid phase. We then computed 
various many-body correlators in the normal phase 
close to criticality. We showed that deviations from the 
Hartree-Fock theory are important, and they show up 
in the renormalized density-density correlator and the 
reduction of the three-body recombination rate. We 
flnally extended the results to the inhomogeneous case 
of trapped Bose gases. The density proflles for various 
phase-space densities were evaluated and we found that 
close to criticality, the RG approach becomes necessary 
for a quantitative description of the gas. We hope that 
these beyond mean-fleld effects can also be observed 
experimentally in the near future. 
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VI. CONCLUSION 



To conclude, we have studied various aspects of the 
2D ultracold Bose gas. First, we derived the exact 
form of the T-matrix for the 2D system as realized in 
experiments. The 2D effective interaction assumes a 
form which interpolates between the 2D and 3D results. 



APPENDIX A 



In this appendix we describe in detail the derivation of 
Eq. (10). Starting with Eq. (9), we add and subtract the 
quantity J2n I G^^k|(/)n(0)P/(27r)^£^o(k, n) to obtain 



m 



m 



T{E) 47rn?a AttH^ 



n=0 



4n + 1 



4n - 2E/hu;z 



n=0 



(27r)2£;o(k,n) J (27r)3 2eq' 



(Al) 



where the zero-point energy huJz/2 has been subtracted 
from the energy argument and the integration 
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(2^)2 



1 



1 



E + - Eo(2n, k) ^o(2n, k) 

I 



m 



In 



4n+ 1 



47r;i2 L4n-2E/;kj;,_ 



has been carried out. The third term in the right-hand 
side of Eq. (Al) can be written as 



J 



(2^)2^o(k,n) .'^^^J (2^f 



27r/22^(2n)! 



^o(k,2n) 



d^k 



hm , ^ 



g-.V4^ g (-l)-i/2n(^/V2/) 1 



= hm / 



n=0 



2n + z/ 



e--V4^ gL^^/')(xV2/2) 



n=0 



2n + i/ 



r 



(A2) 



where the variable = 1/2 + /^ k^ is conveniently de- 
fined, L^n^^'^\x) is the generalized Laguerre function, 
and the relations iJ2n(0) = (— l)"^(2n)!/n! and H2n{x) = 
(-l)^22^n!Li"^/^^(x2) have been used. 

Following Ref. [31], we use the integral representation 



2n 



1 r 

Jo 



dy 



y 



2n+v>-l 



for 2n + > and one of the generating functions of the 
generalized Laguerre function 

f:4-^/2)(x>" = (l-.)-V2exp(-^). 



n=0 

In this manner we find 



Jo 



dy 



v-2 



T{l-v/2YV{v) sin2(7rz//2) 



The last two terms in the right-hand side of Eq. (Al) can 
then be numerically integrated to give 



oo ^ 



(27r)2i;o{k,n) J (27r)3 2eq 



/ 



d^ci 1 



r3/2 



m f 



r d^k 

J (27r)2r(l 



2^ 



(l-z//2)^r(z/) sin^(7rz//2) 



d^k 1 
(2^)2 |k| 



-0.0279 m/Z/i^ 



(A3) 



Finally, by using roo(£^) 
tained. 



3(0)pT(£^), Eq. (10) is ob- 
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